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1.  Introduction 


The  renormalization  group  (RNG)  procedure  for  nonlinear,  dissipative  sys¬ 
tems,  introduced  by  Wilson  (1975),  has  since  been  applied  to  the  Navier  Stokes 
equations  (Forster,  Nelson  and  Stephen,  1978,  De  Dominicis  and  Martin,  1979, 
Fournier  and  Frisch,  1983,  and  Yakhot  and  Orszag,  1986).  Dynamical  equations 
for  the  large  scales  and  long  times  are  derived  by  averaging  over  an  infinitesimal 
band  of  small  scales  to  remove  the  small  scales  from  explicit  consideration.  This 
procedure  yields  infinitesimal  modifications  in  the  equations  for  the  large  scales. 
The  removal  process  is  iterated  and  these  corrections  accumulate  to  give  finite 
changes.  In  the  case  of  the  Navier  Stokes  equations,  the  resultant  model  equa¬ 
tions  for  the  large  scales  are  characterized  by  an  eddy  viscosity  that  replaces  the 
molecular  viscosity. 

The  present  formulation  of  the  RNG  theory  assumes  that  inertial-range  tur¬ 
bulence  can  be  modeled  by  isotropic  turbulence  driven  by  a  random  force.  The 
force  represents  the  net  effect  of  the  energy  cascade  from  the  very  large  scales  on 
eddies  in  the  inertial  range,  and  is  chosen  to  reproduce  the  Kolmogorov  energy 
spectrum  when  opposed  by  the  modified  viscosity.  Homogeneity  in  space  and  time 
and  isotropy  in  space  allow  transformation  to  Fourier  space,  where  the  equations 
are  algebraic. 

For  the  purpose  of  large-eddy  simulations,  one  would  like  to  derive  an  eddy 
viscosity  for  all  wavenumbers  0  <  fc  <  A,  where  A  is  the  cutoff  wavenumber 
defining  the  resolution  of  the  simulation.  Assuming  a  purely  eddy  viscosity  model, 
other  studies  of  isotropic  turbulence  (Kraichnan,  1976,  Chollet  and  Lesieur,  1981, 
and  Domaradzki,  Metcalfe,  Rogallo  and  Riley,  1987)  indicate  that  a  realistic  eddy 
viscosity  exhibits  a  plateau  value  far  from  the  cutoff  (k  «  A)  and  increases  to 
several  times  its  plateau  value  very  near  the  cutoff  {k  «  A).  This  “cusp-up” 
behavior  of  the  eddy  viscosity  is  predicted  by  the  Test-Field  model  of  Kraichnan 
(1976)  and  EDQNM  (Chollet  et  al.  1981),  and  is  consistent  with  direct  numerical 
simulations  (Domaradzki  et  al.  1987). 

At  each  iteration  of  the  RNG  fine-scale  removal,  the  correction  to  the  viscosity 
is  found  in  the  low- wavenumber  limit  k  —*  0.  Thus  the  RNG  subgrid-scale  model 
(Yakhot  and  Orszag,  1986)  includes  the  traditional  eddy-viscosity  type  triads. 
These  triads  have  k  «  A,  and  p,q  >  A,  where  k  -f  p  -1-  q  =  0.  As  discussed  by 
Kraichnan  (1976),  the  triads  leading  to  the  cusp  are  of  a  different  nature,  with 
k  just  below  the  cutoff  A,  p  <<  A  and  q  just  above  A.  This  work  explores  the 
possibility  of  relaxing  the  limit  k  —*  O  in  the  derivation  of  the  RNG  eddy  viscosity. 

In  fact,  the  RNG  procedure  involves  the  two  limits  k  0  and  SA  — >  0,  where 
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SA  is  the  width  of  the  removed  high-wavenumber  shell,  and  we  have  determined 
that  these  two  limiting  processes  do  not  commute.  The  eddy  viscosity  derived  by 
Yakhot  and  Orszag  (1986)  is  found  for  the  ratio  k/SA  — »  0,  and  we  show  that 
k/SA  — »  0  is  necessary  to  recover  physically  meaningful  results. 

Since  SA  is  an  unphysical  parameter,  we  believe  that  the  restriction  k/SA  —*  0 
must  be  re-interpreted  as  a  condition  on  a  Reynolds  number,  and  may  be  equiv¬ 
alent  to  the  so-called  e-expansion.  It  is  possible  that  the  cusp-up  can  be  derived 
from  extrapolation  of  an  expansion  of  the  effective  viscosity  near  k  0  (for  exam¬ 
ple  a  Fade  scheme).  However,  it  seems  likely  that  the  interactions  that  lead  to  the 
cusp  behavior  in  a  purely  eddy  viscosity  model  are  represented  by  different  terms 
in  the  RNG  model,  for  example  Galilean  invariant  higher-order  nonlinearities. 

Also  while  at  Universite  Libre  de  Bruxelles,  the  principle  investigator  under¬ 
took  to  derive  the  vorticity  equation  for  long  times  and  large  scales  using  the  RNG 
method.  The  results  are  described  in  the  enclosed  preprint  which  will  be  submit¬ 
ted  as  a  brief  communication  to  Physics  of  Fluids  A.  The  restriction  k/SA  0  is 
implicitly  enforced,  and  the  long-time,  large-scale  vorticity  equation  is  compared 
to  the  long-time,  large-scale  velocity  equation.  Information  about  the  large-scale, 
long-time  vorticity  equation  should  be  helpful  for  large-eddy  simulations. 

The  second  section  of  this  report  reviews  the  RNG  procedure  to  remove  small 
scales.  The  third  section  shows  that  iterative  removal  of  scales  from  the  Navier- 
Stokes  quadratic  nonlinearity  gives  the  Yakhot-Orzsag  eddy  viscosity  for  k/SA  — ^ 
0,  and  zero  for  SA/k  — »  0.  In  the  fourth  section,  we  attempt  to  clarify  the  work 
of  Zhou,  Vahala  and  Hossain  (1988,  1989).  In  an  attempt  to  derive  the  cusp, 
they  considered  iterative  removal  of  scales  from  the  cubic  nonlinearity  that  arises 
after  removal  of  one  shell  from  the  original  quadratic  nonlinearity.  We  show  that 
the  results  for  k/SA  — »  0  are  not  mathematically  meaningful.  Carati  (1991) 
previously  showed  that  the  results  are  unphysical  for  SA/k  —>■  0,  and  we  give  a 
review  of  Carati’s  (1991)  work.  In  Section  5,  it  is  suggested  that  the  source  of 
problems  with  the  cubic  term  is  non-Galilean  invariance  in  the  limit  of  long  times. 
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2.  The  RNG  scale- removal  procedure 

The  RNG  model  is  for  a  Newtonian  fluid  in  an  infinite  domain,  stirred  by  a 
Gaussian  random  force  f , 


+  VjVjVi  =  -Vip  +  +  fi 


ViVi  =  0  (2.1) 

where  Vo  is  the  molecular  viscosity  and  the  density  p  has  been  absorbed  into  the 
pressure  p.  The  force  is  defined  by  its  two-point  correlation  function  in  wavevector 
and  frequency  space, 


<  /.•(k,a;)/,(k',  w')  >=  2Do{2T^Y^^k-yPij{\!i)S{\i  -f  k')6(u;  -f  w'),  Al  <  k  <  Ao 


=  0,  0  <  A;  <  Ai, 


(2.2) 


where  k  is  the  wavevector,  uj  is  the  frequency,  k  =  [k]  and  d  is  the  number  of 
space  dimensions.  The  parameter  y  is  chosen  to  give  the  Kolmogorov  form  of  the 
energy  spectrum  and  in  three  dimensions  y  =  d  =  3  (Yakhot  and  Orszag,  1986). 
The  wavenumber  Ao  is  aji  ultraviolet  cutoff,  beyond  which  the  energy  spectrum 
is  taken  to  be  zero  and  the  viscosity  is  the  molecular  viscosity  i/g.  Statistical 
homogeneity  in  space  and  time  is  guaranteed  through  ^(k  -|-  k')6(u;  -t-  u;').  The 
projection  operator  Pij(k.)  =  Sij  —  kikjfP  makes  the  force  statistically  isotropic 
and  divergence-free.  In  the  limit  of  an  infinite  Reynolds  number,  the  integral  scale 
Al->0. 


Small-scale  fluctutations  are  eliminated  from  explicit  appearance  in  (2.1)  by 
averaging  over  the  force  field  at  small  scales.  The  force  given  by  (2.2)  is  postulated 
to  represent  the  end  result  of  the  energy  cascade  from  large  to  small  scales,  and 
therefore  models  the  production  due  to  the  quadratic  nonlinearity. 


The  Fourier  transform  of  (2.1)  leads  to 


«i(k)  =  G(k)/i(k)  -  5^G(k)P™„(k)  J  ti„(q)0„(k  -  q)  (2^  (2-3) 
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where  <  fc  <  A^,,  k  =  (k,a;),  G(k)  =  (-iu>  +  ^imn(k)  =  fc„Pm(k)  + 

^n-ftm(k)  and  dq  denotes  the  (d  +  l)-dimensional  integral  over  the  wavevector 
and  frequency  components  of  q.  The  parameter  =  1  is  simply  an  ordering 
parameter  for  the  dimensional  equations,  used  to  identify  where  the  Reynolds 
number  would  appear  in  the  nondimensional  form  of  (2.3). 

•  To  begin  the  RNG  scale- removal,  one  first  defines  the  high- wavenumber 
band  below  A^,  to  be  removed,  A  <  k  <  A,,,  where  the  shell  width  is  denoted 
SA, 

M  =  (A,A,].  (2.4) 

Then,  one  separates  the  velocity  modes  into  v'^  and  according  to  scale, 

v<(k)  =  v(k),  k<A 

.  =  0,  A<k<Ao  .  . 

v>(k)  =  0,  k<A  ^  ^ 

=  v(k),  A  <  k  <  Ao. 

Then  equation  (2.3)  is  written 

«,(k)  =  G(k)/.(k)  -  5^G(k)/>i„»(k)  I (k  -  q)  + 

2ii>(q)«„<(k  -  q)  +  t;>(q)6„>(k  -  q))^^.  (2-6) 

•  Since  the  local  Reynolds  number  of  the  modes  in  the  shell  is  small,  one 
expands  in  a  perturbation  series  in  powers  of  Xg 

v>  =  vW  +  A<,v<*)  -I-  A2v(2)  0{Xl) 

One  finds 

i;?(k)  =  G(k)/?(k)  +  ^G(k)i'™„(k)  J (v<(q)0i(k  -  q)  + 

2G{q)/>(q){i;(k  -  q)  +  G(q)/>(q)G(k  -  q)/>(k  -  «)  (2-7) 

keeping  in  mind  that  A,,  identifies  where  the  local  Reynolds  number  appears 
in  the  nondimensional  equations.  The  velocity  scale  for  the  local  Reynolds 
number  is  the  square  root  of  the  energy  beyond  the  cutoff  Ag,  and  the  time 
scale  is  based  on  the  viscosity  at  the  cutoff. 
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•  One  substitutes  the  perturbation  expansion  for  v^(q)  into  the  equations 
for  the  velocity  modes  v<(k)  at  large  scales  k  ^  Az,  0  and  long  times 
a;  — >  0.  An  ensemble  average  over  the  force  field  in  the  high- wavenumber 
band  using  (2.2)  eliminates  the  small  scales  from  the  equations  for  the  large 
scales.  The  large-scale  variables  are  assumed  independent  of  the  ensemble 
average  over  the  force  at  high  wavenumbers,  and  this  condition  is  nearly 
satisfied  for  k  — ^  0.  The  resulting  equation  for  are  written  symboliccilly, 
without  subscripts,  arguments,  or  coefficients,  as 

(-iu:+h(k,A))v<  =  f<+F<+Xo  J(v<f+Xl  j{v^f+0{Xl  j  {v<Y).  (2.8) 
where  h{k^Ao)  =  Vok^-, 


h{k,A)-h{k,Ao) 


XlDo  Pib(k) 


Plmnik)  J 


Pnabjk-  q)Prr,aiq)q  ^dq  4 

(2.9) 


and  F'^  is  an  induced  force  acting  at  low  wavenumbers  (Yakhot  and  Orszag, 
1986).  The  frequency  integration  over  fl  has  already  been  carried  out  using 
the  Residue  theorem,  where  q  =  (q,  0),  and  we  have  taken  the  limit  w  — >  0, 
where  k  =  (k,a;).  The  wavevector  integration  in  (2.9)  is  restricted  by  A  < 
q  <  Ao  and  A  <  |k  —  q|  <  Aq.  The  origin  of  the  0(Xl)  term  in  (2.9)  is 
the  r)^(q)u^(k  —  q)  term  in  (2.6),  and  the  implications  of  this  are  discussed 
below.  Yakhot  and  Orszag  (1986)  took  the  limit  kfq  0  (also  implicitly 
assuming  k/SA  — >  0,  see  Section  4)  and  found 


hik,A)  -  h(k,A,)  =  MA)  -  UT{A,))k'^ 


t'r(A)  —  i'r(Ao)  = 


AAl  DoSd 

1/2  (27r)‘^  (4  -f  y  -  d) 


(2.10) 


(2.11) 

where  Sd  is  the  area  of  a  sphere  in  d  dimensions.  The  reason  for  the  intro¬ 
duction  of  the  new  parameter  e  will  become  clear  below. 


Iteration  of  the  scale  removal  procedure  can  be  performed  by  integration  of 
the  differential  equation  for  vt{A), 
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(2.12) 


..  J^(A)  —  i/3’(Ao)  _  Adt^{A)X^{A) 

dA  ^  «A-o  Ia  “  “”2  A 


+  O(A^) 


P(A)  = 


2DoSd  XI 

{2tY  4(A)A‘’ 


(2.13) 


The  quantity  A(A)  is  the  effective  Reynolds  number  at  A,  based  on  the  modified 
viscosity  vt{A). 

Notice  that  in  writing  the  differential  equation  (2.13)  one  has  assumed  that 
the  value  of  uj  derived  in  the  limit  fc  — »  0  is  valid  for  all  k  up  to  the  cutoff  A 
(Smith  and  Reynolds,  1991,  Section  3).  Analysis  of  direct  numerical  simulation 
data  (She,  1991)  shows  that  the  effective  eddy  viscosity  due  to  triad  interactions 
between  wavevectors  k,  p  and  q,  where  k  +  q  +  p  =  0,  k<A  and  p,  q  >  A,  is 
indeed  flat  for  0  <  A:  <  A  (She,  1991).  Thus  this  assumption  is  consistent  with 
the  origin  of  the  RNG  eddy  viscosity  being  the  ti>(q)t)>(k  —  q)  term  in  (2.6). 

With  initial  condition  t'r(Ao)  =  t'o?  integration  of  (2.12)  leads  to 


vt{A)  = 


ZAdD,Sd{A-^-A-^)^^l^ 

{2irY  c  J 


(2.14) 


with  (2.14)  strictly  valid  for  A  <<  Ao  where  it  is  independent  of  Vo  and  Ao  in 
keeping  with  the  assumed  similarity  range.  For  A  <<  Ao, 


(2.15) 


and  by  (2.13) 


For  e  -+  0  and  using  (2.11),  one  sees  that 


(2.16) 


+ 


(2.17) 


where  A®  =  (d^  —  d)f{2d{d  +  2)).  Therefore,  for  e  — ^  0,  (2.12)  becomes 


(2.18) 


dvT{\)  (^3  +  0(e))i/T(A)A^(A) 


dK 


+  0(cA2). 


For  a  consistent  asymptotic  expansion,  all  O(cA^)  terms  must  be  dropped,  and 
the  lowest-order  solution  for  j^(A)  is 


strictly  valid  in  the  limit  A  <<  Ao-  Thus  for  c  0  and  A  <<  A^, 


Iw)  • 


(2.20) 


Since  A  is  the  effective  Reynolds  number  of  the  flow  at  A,  all  perturbation 
expansions  used  in  the  scale-removal  procedure  are  exact  for  c  — v  0.  Furthermore, 
the  modified  equation  for  v'^(k)  is  given  at  lowest-order  in  cin  expansion  in  powers 
of  c.  The  Yakhot-Orszag  theory  is  an  extrapolation  of  the  RNG  results  at  lowest 
order  in  t  to  the  case  c  =  4,  which  is  the  value  for  Kolmogorov  scaling  (Yakhot 
and  Orszag,  1986). 

Notice  that  the  cubic  nonlinearity  in  the  nondimensional  form  of  (2.8)  is 
O(A^)  =  0(e)  and  can  be  neglected  in  comparison  to  the  quadratic  nonlinear¬ 
ity,  which  is  0(A)  =  0(e^^^).  However,  Zhou  et  al.  (1988,  1989)  considered  the 
contribution  to  h{k,  A')  after  removal  of  the  next  band  A'  <  A:  <  A,  and  fur¬ 
thermore  the  accumulated  contribution  from  the  cubic  after  iterative  removal  of 
scales.  Their  goal  Wcis  to  derive  the  cusp  in  the  eddy  viscosity  for  A:  «  A,  where  A 
is  then  the  final  cutoff,  since  the  cubic  term  arises  from  the  t)^(q)t)^(k  —  q)  term 
in  (2.6).  Kraichnan  (1976)  showed  that  the  cusp  is  a  result  of  triad  interactions 
with  k  just  below  the  cutoff,  q  just  above  the  cutoff,  and  p  much  smaller  than 
the  cutoff.  This  has  been  verified  by  analysis  of  direct  numerical  simulation  data 
(She,  1991). 

In  Section  4,  we  show  that  the  contribution  to  h(k,  A')  does  not  make  mathe¬ 
matical  sense  for  k/SA  —>■  0,  which  condition  is  necessary  to  derive  the  eddy  vis¬ 
cosity  from  the  quadratic  nonlinearity,  and  is  probably  related  to  the  c-expansion. 
Carati  (1991)  showed  that  the  contribution  to  h(k,A')  from  the  cubic  is  unphys¬ 
ical  in  the  case  6A/k  —*  0.  A  fundamental  problem  with  the  cubic  is  that  it  is 
not  Galilean  invariant  in  the  limit  of  long  times  (a;  — »  0),  and  this  is  probably  the 
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reason  why  the  equations  including  the  cubic  do  not  give  meaningful  results  when 
8A/k-*0. 

Section  3  shows  tl^at  the  condition  k/8A  — +  0  is  necessary  to  derive  the  RNG 
eddy  viscosity  from  the  Navier  Stokes  quadratic  nonlinearity.  If  this  condition 
is  relaxed,  there  is  no  contribution  to  h(k,A)  from  the  quadratic  nonlinearity. 
Section  4  shows  that  the  cubic  nonlinearity  yields  nonsensical  results  for  k/8A  — » 
0.  It  yields  unphysical  results  for  8A/k  — »  0,  probably  because  it  is  not  Galilean 
variant  in  the  limit  of  long  times  (Section  5). 
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3.  The  renormalized  viscosity 

The  equation  for  v'^  is,  from  (2.6) 


G-HHk)  =  /i'(k)  +  Jlv<(q)  v„<(k  -  q)+ 

2<>S(q)i’J(k-q)  +  tiJ(q)«?(k-q)l  <*■') 

where  G~^  =  {—iiv  +  h{k,\)).  Substituting  the  A-exp^insion  (2.7)  for  v^,  up  to 
second  order  in  Aq,  generates  two  terms  linear  in  v*'  and  quadratic  in  /^.  One 
comes  from  the  v^v'^  term  but  it  vanishes  after  ensemble  averaging  over  and 
integration  over  q.  The  other  linear  term  is  generated  from  the  term  in  (3.1), 
symbolically  this  contribution  is  2(G^/^)(Ao  2G^  Ensemble  averaging  over 

/>  using  (2.2),  the  explicit  form  of  this  second  linear  term  is,  from  (3.1)  and  (2.7), 

2Dq\1  J  P,m„(k)Pmr,(q)/'nr(k-q)  X 

G(k-q)G(q-k)G(q)  |k  -  «.(k)  (3.2) 

where  the  integration  domain  is  such  that  A'  <  |k  —  q|,  q  <  A.  This  is  the 
expression  (2.9).  This  expression  has  the  form  7i>(k)u,(k),  but  from  isotropy  the 
tensor  Tu(k)  must  be  proportional  to  P,-,:  T»»(k)  =  6G~^  -ft»(k)  where  the  scalar 
6G~^  =  Pia{k)Ti3{k) / {d  —  1).  Hence  this  contribution  can  be  combined  with  the 
propagator  G(k,  A)  to  define  a  new  propagator: 

{G')-\kA')  =  G-\k,A)  +  6G-'^ 


where 

2M  r 

SG-^  =  (^  /  P>mnik)Pmrs{q)Pnr{p)  G(p)G(-p)G(q)  p-^q  (3.3) 

with  p  =  k  —  q  and  A'  <  p,  q  <  A.  The  factor  M  stands  for  M  =  Dq/KH  — 
l){2nY-^]. 

Taking  the  long-time  limit,  integrating  over  Izirge  frequencies  fl  (where  q  = 
(q,  fi))  using  the  Residue  theorem,  and  making  the  usual  change  of  variable 

jdq  =  jy<l,j  j?ldpd. 
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where  <f>  determines  the  orientation  of  the  triad  around  k,  the  equation  for  h{k,  A) 
becomes,  from  (3.3) 


h{k,A')-hik,A)  =  M  j 

J  D‘. 


A{k,p,q)p  " 


/m  h{p.  A)  [A(p,  A)  +  h{q,  A)]  k 
The  integration  domain  D2  is  such  that  A'  <  p,  9  <  A,  and 

Mk,P,q)  =  i’»mn(k)P„r,(q)Pmr(p)  =  “  P^)(9^  -  P^)  + 

with  =  2k^p^  +  2k^q^  +  2p^9^  —  k*  —  p*  —  q*. 


^dpdq  (3.4) 


Ak'^p'^q'^ 


3.1.  k  <  6A 

When  k  <  SA,  the  domain  of  integration  for  the  contribution  from  the  quadratic 
term  (darkest  area  in  Fig.  1),  is  composed  of  three  distinct  parts  in  the  p,  ^-space, 


r  /A  fS.  f\-k  tP+k  r 

/  =  /  dp  dq+  dp  dq  + 

JD2  J\-k  JA-k  Ja'  Jp  JA 


fA—k 

/A' 


Fig.  1.  Domain  of  integration  when  k  <  6 A. 
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When  A  =  Ao,  h{k,  A)  =  vq  where  i/q  is  the  molecular  viscosity.  Assume 
that  A)  =  i/{A.)k^,  for  all  k.  Expanding  in  powers  of  k  (with  9  =  p  + 
kx,  — 1  <  ar  <  1)  and  6A  =  A— A',  one  finds  after  some  lengthy  but  straightforward 
manipulations, 


M 

h{k,  A')  -  h(k,  A)  =  -k^  8v  =  X 

i/^(A) 

[  — ib*  8k  +  P  A-^-"  +  k-^-y  8k  +  0{k\8k^)] 

In  the  limit  k/8k  0,  one  obtains  a  differential  equation  for  t'(k)  =  h{k,k)/k'^ 

which  is  readily  integrated  to  find 

.-“(A)  -  ,/“(A„)  10  -  2y  A-'-»  -  A„-'-» 

— 3 —  = "  ^5 - (Tin;)— 

This  expression  is  identical  to  that  (2.14)  obtained  by  Yakhot  and  Orszag  (1986) 
when  d  =  3  and  y  =  e  —  \  (only  d  =  3  is  considered  in  this  report). 

3.2.  k>  8k 


Figure  2.  shows  clearly  that  when  k  >  8k  the  contribution  from  the  quadratic 
term  is  second  order  in  ^A.  The  domain  of  integration  D2  (darkest  area  in  Fig.2) 
in  (3.4)  is  simply 

/  =fdpfdq 
JD2  Ja'  Ja’ 

then  (3.4)  with  h{k,k)  =  i/(A)P  gives 


h{k,k-8k)-h{k,k)  = 


M  A{k,k,k)  k-y  k? 

2 


4A‘*  k 
M  A-'‘-5'(4Ait  -  P) 
1/2  8 


^(6A)2  +  0(<5A)3 
{8kf  +  0{8kf 


One  obtains 

k^  8u  oak  (dA)^ 

hence  the  quadratic  term  does  not  contribute  in  the  limit  6A  — »  0  for  fixed  k. 
Note  that  the  form  of  the  increment  (oc  k)  is  not  consistent  with  the  assumed  k^ 
dependency  for  h{k,k). 
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Fig.  2.  Domain  of  integration  when  k  >  SA. 


4.  The  cusp  and  the  triple  non-linear  terms 

The  subgrid  scale  eddy- viscosity,  i/{k/A),  sis  defined  by  Kraichnan  (1976), 
shows  a  very  steep  rise  (a  “cusp”)  as  k  approaches  the  numerical  cut-off  A.  As 
explained  by  Kraichnan,  this  cusp  behavior  is  due  to  the  effect  on  a  wavenumber 
just  below  the  cut-off  of  triad  interactions  with  a  very  small  wavenumber  and 
a  wavenumber  just  above  the  cut-off.  In  the  language  of  the  previous  sections, 
this  would  originate  from  the  2v'^u^-term  in  the  -equation.  The  sharp  cusp  in 
the  eddy-viscosity  shows  that  this  effect  is  not  well-modeled  by  an  eddy-viscosity. 
Note  also  that  the  cusp  is  not  universal  but  depends  on  the  large  scales  (Kraichnan 
1976). 

After  one  step  of  the  shell-elimination  procedure,  the  2v'^u^-term  generates  a 
triple  non-linear  term  which  is  not  present  in  the  original  Navier-Stokes  equations. 
(It  is  shown  below  that  the  triple  term  is  not  Galilean  invariant).  However,  at  the 
next  shell-elimination  this  triple  non-linear  term  generates  a  new  linear  term.  Zhou 
et  al.  (1988,1989)  proposed  that  this  new  contribution  to  the  propagator  could 
be  the  one  responsible  for  the  cusp  behavior  of  the  subgrid  scale  eddy-viscosity. 
They  went  on  to  implement  a  numerical  version  of  the  iterative  shell-elimination 


procedure.  Unfortunately,  their  results  are  strongly  dependent  on  the  shell-width 
(6A)  selected  for  the  numerics.  There  is  no  evidence  that  the  results  converge 
when  the  shell- width  tends  to  zero.  To  clarify  this  issue,  Carati  (1991)  studied 
the  problem  in  the  limit  of  infinitesimal  shells.  He  showed  that  the  contribution 
from  the  triple  term  to  the  propagator  is  proportional  to  k  for  small  k  and  not 
to  k^  as  for  a  viscous  term.  This  k  contribution  does  not  seem  to  make  physical 
sense. 

The  shell-elimination  procedure  (Sect.2)  starts  by  introducing  a  cutoff  Ai  < 
Ao,  and  splitting  the  velocity  field  v  into  v'^  and  v^.  Substituting  the  A-expansion 
(2.7)  for  nm(Q)  in  the  u^u'^-term  of  equation  (3.1),  generates  a  triple  non-linearity 
on  the  right-hand  side  and  a  term  linear  in  (Sect. 3).  The  triple  non-linearity 
has  the  form: 

j  J  2  -^Go(q)Pmr.(q)  [n,'^(q)  -  Q)  Vn  (k  -  q)] 

(4.1) 

At  the  next  step,  a  new  cutoff  A2  <  Ai  is  introduced  and  the  velocity  field 
(formerly  u'^)  is  separated  once  more  into  <  and  >  components.  One  uses  the  A- 
expansion  for  the  new  terms.  The  triple  non-linearity  generates  various  terms, 
in  particular  a  linear  term  of  the  form 

(  G.(q)G.(q  -  Q)  /’(Q)  /i-(q  -  Q)  c; (k  -  q)  + 
G,(q  -  q)G.(k  -  Q)  «?(q)/,>(q  -  Q)A>(k  -  q)  + 
G,(q)G,(k  -  q)  />(Q)«f(q-  Q)/„=-(k  -  q)  J 

Averaging  over  realizations  of  the  force  with  (2.2) 

<  Mk)fj{k')  >=  2Do{2Ty+^k-y6{k  k')^.i(k) 

the  term  proportional  to  Pmrs(q)  <  />(Q)/>(q-Q)  >  vanishes  upon  integration 
over  Q,  q.  One  is  left  with 

2B.  Xl  J  «„„(k)/>„,.(q)F„.(p)  G,(p)G,{-p)G„(q) 

where  p  =  k— q.  This  expression  is  quite  similar  to  that  obtained  from  the  original 
quadratic  non-linearity  (3.2)  and  one  proceeds  as  before.  Taking  the  long-time 
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limit  (a;  — +  0)  and  integrating  over  large  frequencies,  one  defines  a  new  propagator 
at  A2: 


Do>^l  f  /^.mn(k)Pmr»(q)Pnr(p)p~" 

{d-l){2iryJ  fi(p,Ai)[/i(p,Ai)  + /i(g,Ao)]  ^ 


The  differences  with  expression  (3.4)  reside  in  the  domain  of  integration,  which 
is  not  symmetric  in  p,  q,  and  the  A-dependency  of  the  functions  h{k.  A)  in  the 
integrand. 

Of  course  a  new  triple  non-linearity  has  been  generated  after  this  second 
shell-elimination  and  it  will  also  contribute  to  the  propagator  at  the  next  shell- 
elimination.  After  A^  -f  1  eliminations,  the  sum  of  all  contributions  from  the  triple 
non-linearities  to  the  propagator  can  be  written: 


h(k,  —  h(k,  An)  =  M  VI  f  t7 — l — \  — TT — (4-3) 

^  ^  ^  ^qJb^  h{p,AN)[h{p,AN)  +  h{q,Aj)]  k  ’ 


The  integration  domain  Bf  is  such  that  An+x  ^  P  ^  and  Aj+i  <  q  <  Aj. 


4.1.  k<8A 

Suppose  that  all  wavenumbers  down  to  A  have  already  been  eliminated  and  con¬ 
sider  now  the  next  step  of  the  procedure  where  another  shell  of  wavenumbers  will 
be  eliminated  to  reach  A'  <  A.  When  k  <  SA  there  is  only  one  term  in  the  sum 
(4.3),  that  where  An+i  <  P  <  An  and  An  <  9  <  Aat-i.  The  other  domains 
are  empty  (Fig.l).  The  domain  of  integration  D3  =  (lightly  shaded  area  in 
Fig.l)  for  the  cubic  term  is  independent  of  A'  when  k  <  6A,  it  is  given  by 


f  fA+k  rA 

I  =  I  dq  I  dp 
JD3  Ja  Jq-k 


Expanding  in  powers  of  k  and  SA,  as  before,  the  contribution  of  the  cubic  to 
h(k,  A')  —  h{k,  A)  =  —k^S(c)^  is  equal  to 


(4.4) 


This  contribution  does  not  make  mathematical  sense  as  the  right-hand  side 
should  be  of  order  k'^SA  but  is  of  order  k"^  only. 
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4.2.  k>6A 


This  case  has  been  studied  in  detail  by  Carati  (1991).  The  analysis  is  sketched 
below.  In  the  limit  of  infinitesimal  shells,  Aj+i  —  Aj  =  6Aj  — 0,  N  oo,  letting 
Ajv  =  A,  the  equation  (3.8)  becomes 

aHk,\)  _  MA'-»  /■*«  A,A(t,A,Aj) 

dA  kh(A,A)  h  |A(A,A)  +  A(A,,Aj)]  •'  '  ' 

This  is  an  integro-differential  equation  for  h{k^A).  Carati  (1991)  looked  for  a 
self-similar  solution  of  the  form 


/i(jb,A)  =  AV(A:/A)  (4.6) 

Substituting  in  (4),  one  finds  c  =  (5  —  j/)/3  for  self-similarity  and 

with 

F{v)  =  v  J 


To  determine  the  form  of  the  propagator  (4.6)  in  the  limit  of  very  large  scales, 
one  has  to  expand  <f{l)  around  /  =  0.  Let 

OO 

F(v)  =  -£F>„< 

i=0 


then  integrating  (4.7)  term  by  term  with  v  =  lx  yields 


From  (4.8),  one  finds 


V  .  4  -t-  c  2  .  26  —  3c 


Writing  c  =  (5  —  j/)/3  =  2  —  e/3,  the  final  result  reads 
'■(*^,A)  =  ^A-/^x 


.  -3  k  3  18  -  e  P  3  20  +  t  P  A: , 


(4.9) 


(4.10) 


(4.11) 


The  lowest  order  term  is  proportional  to  k  and  this  term  does  not  have  any  clear 
physical  interpretation.  It  does  no'  .  .  re'-nond  to  a  renormalized  viscous  damping. 
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5.  Non-Galilean  invariance  of  triple  non-linearities 

In  this  section  it  is  shown  that  the  triple  non-linear  term  generated  by  the 
shell-elimination  procedure  is  not  Galilean  invariant.  It  is  believed  that  this  non¬ 
invariance  of  the  resulting  equations  is  the  cause  for  the  non-sensical  results  pre¬ 
sented  in  Sect. 4.  The  reason  is  that  at  the  second  cut  one  performs  the  shell- 
elimination  on  an  equation  which  has  lost  an  essential  symmetry  of  the  original 
Navier-Stokes  equations,  the  Galilean  invariance.  Thus  keeping  the  triple  non- 
lineajity  is  in  conflict  with  the  spirit  of  the  Renormalization  Group  analysis  of  the 
Navier-Stokes.  At  each  step  of  the  RNG  procedure,  the  equations  are  renormal¬ 
ized  to  obtain  new  equations  formerly  identical  to  the  original  equations,  and  in 
particular  preserving  all  the  symmetries  of  those  original  equations. 

Consider  a  Galilean  transformation  : 

X  =  x'  -h  V  t 
t  =  t' 

v(x,t)  =  V-hv'(x',t') 

In  Fourier  space: 

k  =  k' 
a;  =  u'  +  k'-V 
v(k,a>)  =  V(5(k)^(a;) -f  v'(k,a;') 

As  is  well-known,  a  Galilean  transformation  of  the  Navier-Stokes  equations  (2.3) 
leaves  them  unchanged.  Omitting  viscosity  and  the  force,  which  are  irrelevant  for 
our  purpose  here  (the  white-noise  nature  of  the  force  assures  Galilean  invariance 
on  average),  (2.3)  reads 

-iiO  Vi{k)  =  -^Pimn(k)  j  Vm(q)Vn(k  -  (5-3) 

The  Galilean  transformation  (5.2)  of  (5.3)  leads  to 
^  k'  •  V)  i);(kO  =  -t7’>„„(k')I4v;(kO  + 

/ «„(q')e:(k'  -  q')(5^  (5.4) 

This  has  the  same  form  as  (5.3)  because  Fl>„„(k')y,„v^(k')  =  (k  •  V) 

After  a  first  elimination,  a  triple  non-linearity  of  the  form  (cf.  (4.1)); 

j  J  dqrfQ  P„„(q)  [«?(q)6,‘=(q-Q)«J(k-q)l  (5.5) 


(5.1) 

(5.2) 
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is  added  to  (5.3).  A  Galilean  transformation  (5.1)  of  that  triple  non-linearity  (5.5) 
transforms  it  into,  after  some  reductions  using  the  definition  of  the  function  Pimn 
and  continuity: 


/  j  6i(k'  -  q')  r^q'  -  Q')  + 


i>™„(k')  /  <iQ'  (k-  ■  V)  «;(k  -  Q)  + 
2  fi„„(k)  /rfq  (q  ■  V)  i)i(k  -  q)  {^(q)  + 


2(k  •  v)^e;(k') 


The  appearance  of  the  last  three  terms  makes  the  equations  non-invariant  in  the 
presence  of  the  triple  non-linearity.  The  third  term  is  proportional  to  k. 

Note  finally  from  (5.2)  that  the  limit  of  large  times  a;  — >  0  seems  to  require 
that  the  large  scale  limit  fc  ^  0  be  taken  simultaneously,  in  order  for  the  result  to 
be  Galilean  invariant,  i.e.  a;'  =  a;  —  k  •  V  — »  0,  k'  — >  0,  when  a;  — »  0  and  >  0. 
This  might  suggest  why  the  constraint  kjSti  — »  0  is  required  to  obtain  sensible 
results  using  the  renormalization  procedure. 
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6.  Conclusions 


The  analysis  presented  above  seems  to  indicate  that  the  RNG  procedure  is 
unable  to  derive  the  cusp  in  the  eddy-viscosity,  at  least  at  the  lowest  order  of  the 
technique.  However,  it  should  be  pointed  out  that  this  might  be  more  a  problem 
of  interpretation  thaji  a  shortcoming  of  the  procedure.  The  cusp  arises  when  one 
insist  on  representing  the  effect  of  the  subgrid  sc3des  on  the  resolved  scales  as  an 
eddy-viscosity.  Two  remarks  are  in  order. 

First,  it  is  very  plausible  that  the  RNG  procedure  correctly  determines  that 
part  of  the  interactions  with  the  subgrid  scales  which  can  actually  be  modeled  by 
an  eddy-viscosity.  This  part  would  be  the  only  one  of  importance  when  k  is  less 
than  about  A/2.  When  eliminating  all  wavenumbers  down  to  A  =  fc  while  still 
preserving  the  form  (thus  ignoring  any  cusp),  the  final  result  should  rather 

be  interpreted  as  the  eddy-damping  in  the  inertial  range.  The  RNG  procedure 
can  then  be  used  to  obtain  the  EDQNM  model  as  in  Dannevik  et  al.  (1987)  from 
which  one  can  derive  the  cusp  behavior  in  an  appropriately  defined  eddy-viscosity, 
exactly  as  in  Kraichnan  (1976).  In  effect,  the  RNG  procedure  is  used  in  this  case 
to  determine  the  free  parameter  in  the  EDQNM  model. 

Second,  it  should  be  pointed  out  that  the  cusp  in  the  eddy-viscosity  is  a  spectral 
space  feature.  One  must  include  the  cusp  when  performing  Large  Eddy  Simula¬ 
tions  (LES)  using  spectral  methods  with  a  sharp  cut-off.  It  is  not  clear  at  all  that 
one  needs  to  worry  about  the  cusp  when  using  finite-difference  type  methods.  The 
reason  for  this  is  that  the  cusp  represent  a  statistical  “squishing  out  of  resolution” 
of  small  scales.  The  nature  of  the  conservative  triad  interactions  in  Fourier  space 
makes  it  impossible  for  a  spectral  method  to  let  energy  leak  out  at  the  cutoff  and 
one  observes  an  artificial  pile-up  of  energy  near  the  high- wavenumber  end  if  it  is 
not  removed  by  an  increased  eddy-viscosity  in  that  region.  In  a  finite-difference 
method,  it  is  possible  for  small  eddies  to  be  squished  beyond  resolution,  and  to 
disappear  without  having  to  be  removed  by  extra  amount  of  viscosity. 
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